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Abstract 

We consider constructing the relativistic system of collective coordinates of a 
field theory soliton on the basis of a simple principle: The collective coordinates 
must be introduced into the static solution in such a way that the equation of 
motion of the collective coordinates ensures that of the original field theory. As 
an illustration, we apply this principle to the quantization of spinning motion of 
the Skyrmion by incorporating the leading relativistic correction to the rigid body 
approximation. We calculate the decay constant and various static properties of 
nucleons, and find that the relativistic corrections are in the range of 5% - 20%. 
We also examine how the baryons deform due to the spinning motion. 
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1 Introduction 



Collective coordinates of a field theory soliton play a central role in the study of the soliton 
dynamics (see, e.g., [1] for a review). They represent the motion of the soliton as a whole in 
the symmetry directions of the original field theory. Since the soliton breaks these symmetries, 
the collective coordinates are related with the zero-mode fluctuations around the soliton. The 
simplest way to construct the system of collective coordinates is to promote the (originally 
constant) free parameters of the static soliton associated with symmetries to time-dependent 
dynamical variables by discarding the fluctuation of non-zero modes. 

As an example, let us consider the collective coordinate of space-translation of a static 
soliton (f> c i(x) in a scalar field theory in 1 + 1 dimensions with Lagrangian density C = 
— (l/2)(d fl (f)) 2 — U(4>). In the simple treatment, we take 

0(x,t)=0 cl (x-X(t)), (1.1) 

as the soliton field with a dynamical variable of the center-of-mass X(t). Plugging this into the 
Lagrangian density and carrying out the space integration, we get the following Lagrangian 
of X(t): 

L = -m + imX(t) 2 , (1.2) 

with m being the energy of the static soliton. Somehow, we have obtained the non-relativistic 
Lagrangian of a point particle with coordinate X(t), although we started with a relativistic 
scalar field theory. We need an improved treatment of collective coordinates which leads to 
the relativistic Lagrangian 

L = -mVl-X 2 . (1.3) 

The purpose of this paper is to propose a simple and general principle of constructing the 
system of collective coordinates of a field theory soliton, by which we can obtain the complete 
and relativistic Lagrangian of collective coordinates. 

First, let us state our principle of introducing the collective coordinates into a static soliton 
solution: Collective coordinates must be introduced into the static solution in such a way that 
the equation of motion (EOM) of the collective coordinates ensures that of the original field 
theory. We call this the "EOM principle" hereafter. Although the original field theory has an 
infinite number of degrees of freedom, the reduced system of collective coordinates has only 
a finite number of ones. Therefore, even if the EOM of collective coordinates holds, it does 
not necessarily imply that the field theory EOM holds. Our EOM principle demands in a 
sense that there should be no mismatch between the field theory dynamics and the collective 
coordinate dynamics. If the original field theory is a relativistic one, our EOM principle should 
automatically lead to a relativistic Lagrangian of the collective coordinates. 
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Practically, the introduction of collective coordinates by the EOM principle must be carried 
out in a self-consistent manner. Starting with a suitably chosen dependence of the soliton field 
on the collective coordinates, we obtain the Lagrangian and consequently the EOM of the 
collective coordinates by plugging the soliton field into the field theory Lagrangian. Besides 
this, we plug the assumed soliton field into the field theory EOM, and examine whether it 
holds under the EOM of collective coordinates. We tune the dependence of the soliton field 
on the collective coordinates so that the EOM principle is realized. 

In the above example of the center-of-mass in a scalar field theory in 1 + 1 dimensions, the 
field theory EOM is in fact broken by the term X 2 (p'^(x — X) if we take (1.1) as the soliton 
field (see eq. (1.6)). Let us apply our EOM principle to this system to resolve the breaking 
of field theory EOM successively in the number of time- derivatives. As an improved soliton 
field, we take 4>(x,t) = <f) c i(y(x,t)) with 

y(x, t)=(x- X(t)) (l + c 2 X(t) 2 + c 4 X(t) 4 + •••), (1.4) 

where c 2 and c 4 are constants to be determined by the EOM principle. The corresponding 
Lagrangian of X(t) reads 



/ 



up to time-derivative terms. Therefore, the EOM of X{t) remains unchanged from that in the 
simple treatment: X — 0. On the other hand, for the field theory EOM, we have 

-dl<f>(x,t) + U'(cj>(x,t)) = [(1 - 2c 2 )X 2 + (2c 2 -4- 2c 4 ) X 4 + • • •] ^(y), (1.6) 

up to terms containing X. From this, we find that the EOM principle is satisfied (up to X 6 ) 
by taking c 2 = 1/2 and c 4 = 3/8. Then, y(x,t) (1.4) and L (1.5) are nothing but the first 
three terms of the expansion in powers of X 2 of the Lorentz boost 

/ x x — X(t) „, 
y(x,t) = -=J±L, 1.7 

VI -X 2 

and the relativistic Lagrangian (1.3), respectively* 

Here, we should mention another and conventional way of constructing the system of col- 
lective coordinates. It is, starting with suitably introduced collective coordinates (for example, 
(1.1) for the center-of-mass), to integrate over the non-zero modes around the soliton. Equiv- 
alently, we solve the EOM of the non-zero modes to express them in terms of the collective 



*For a constant X, y(x,t) of (1.7) obviously leads to the relativistic Lagrangian (1.3) to all orders in X 2 . 
However, for a generic X(t), we have to add (1.7) corrections of O(df) starting with |(x — X(t)) 3 (XX) 2 for 
obtaining (1.3). 
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coordinates. Since the set of EOMs of both the zero and non-zero modes are equivalent to 
the field theory EOM, integration over the non-zero modes leads to a system of collective 
coordinates whose EOM implies the field theory EOM. In fact, the relativistic energy of the 
center-of-mass motion, E = \J P 2 + m 2 , has been obtained by this method for a scalar field 
theory in 1 + 1 dimensions [2-4]. Looking back the procedure of (1.4), (1.5) and (1.6), the 
EOM of the non-zero modes implies that they are equal to zero if we adopt the elaborated 
way of introducing the center-of-mass coordinate X(t) by using y(x,t) of (1.4). 

In the case of the center-of-mass coordinate of a soliton, we know its relativistic Lagrangian 
(1.3) from the start. However, for the collective coordinate of the spinning motion of a soli- 
ton, it is quite a non-trivial problem to obtain its relativistic description, and this is our 
main concern of this paper. Concretely, we consider the spinning collective coordinate of the 
Skyrmion, which is a soliton representing a baryon in the low energy effective theory of pions 
(the Skyrme-model) [5]. Let U c \(x) be the static Skyrmion solution. In the simple treatment, 
we take U(x,t) = i7 c i(.R _1 (t)£c) as the spinning Skyrmion field with the rotation matrix R(t) 
of spinning motion. This leads to the Hamiltonian of a spinning spherical rigid body [6,7]: 

frigid body = M cl + ^Ift 2 = M cl + ^/ 2 , (1.8) 

where M cl , X, fl and I(— 2X2) are the energy of the static solution, moment of inertia, angular 
velocity in isospace (body-fixed frame) and the isospin operator, respectively. (H Tigid body also 
has an equivalent expression with fl and J replaced with the angular velocity u) in real space 
and the spin operator J{= Tu>), respectively.) However, this rigid body description is certainly 
a non-relativistic one valid only for a slow spinning motion. For a large angular velocity, the 
Skyrmion should intuitively deform from its original spherical shape, and this should lead 
to corrections with higher powers of f2 to the rigid body Hamiltonian. Therefore, it is an 
interesting theoretical subject to obtain relativistic (i.e., higher time-derivative) corrections 
to the rigid body approximation (1.8), or, if possible, to find the fully relativistic description 
which is the spinning motion counterpart of (1.3). 

Relativistic corrections to the rigid body approximation (1.8) may be important also phe- 
nomenologically for the Skyrmion describing baryons. Identifying the eigenvalues of the 
Hamiltonian (1.8) with the masses of the nucleon (M^ = 939 MeV, I 2 = J 2 = 3/4) and 
A (M A = 1232 MeV, I 2 = J 2 = 15/4), we find the fallowings: 

M cl = 866 MeV, X = 0.00512 MeV" 1 , \Q N \ = 169 MeV, \Q A \ = 378 MeV. (1.9) 

This implies that, (i) about 8% (30%) of the total mass of the nucleon (A) comes from the 
rotational energy, and (ii) the rotational velocity of the nucleon (A) at the radius of 1 fm is 
86% (190% !) of the light velocity. These facts suggest that relativistic corrections are non- 
negligible especially for A, and that we should reexamine the analysis of Refs. [6, 7] based on 
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the rigid body approximation by incorporating the effects of the corrections. The limitations 
of the rigid body approximation were shown numerically without relying on the expansion in 
powers of angular velocity in the case of the (2 + l)-dimensional baby version of the Skyrme 
theory [8]. 

In this paper, we derive the leading relativistic correction to the rigid body Lagrangian 
of the spinning Skyrmion by taking the EOM principle as the basic principle of constructing 
the Lagrangian of spinning collective coordinate. We follow the procedure starting with (1.4) 
for the center-of-mass motion. Namely, assuming the spinning Skyrmion field of the form 
U(x,t) = U c \(y(x,t)) with y(x,t) = + (relativistic correction term), we determine 

the correction term from the EOM principle. Note that the relativistic correction to U (x, t) 
is made only to the argument of the static solution. In the case of the center-of-mass motion, 
the first relativistic correction is specified by the constant coefficient c 2 . For the spinning 
Skyrmion, we have to determine three functions (A(r) , B (r) , C (r)) of the radial coordinate 
r appearing in y(x,t) (see eq. (3.2)) by solving their differential equations derived from the 
EOM principle. Then, from (A(r),B(r),C(r)), the f2 4 correction term to the rigid body 
Lagrangian is obtained (see eq. (4.1)). From the function y(x,t), we can also learn how the 
static Skyrmion of spherical shape deforms due to the spinning motion. 

Then, we repeat the analysis of Refs. [6, 7] for the decay constant f w and various static 
properties of nucleons, such as charge radii, magnetic moments and axial vector coupling, 
by including the relativistic corrections. We find that the contribution of the relativistic 
corrections is in the range of 5% to 20%. However, the comparison with the experimental data 
is a rather disappointing one: Although the value of f n is shifted closer to the experimental 
one due to the relativistic correction, it makes the theoretical value move away from the 
experimental one for most of the static properties. However, this should not be regarded as a 
problem of our EOM principle. Either higher order relativistic corrections are important, or 
we cannot expect the Skyrme model to reproduce the baryon sector so precisely. 

As we mentioned before, there is another way of constructing the system of collective 
coordinates. This is, starting with a simply introduced collective coordinates (such as in 
the rigid body approximation), to solve the EOM of the massive modes to express them 
in terms of the collective coordinates. Though there has appeared no explicit analysis of 
the spinning Skyrmion by this method, the present one using the EOM principle should 
have some advantages: Firstly, ours is simpler since we are not bothered by the non-zero 
modes, and secondly we can directly know how the soliton deforms due to its fast spinning 
motion since the collective coordinates are introduced by deforming the coordinate of the 
static solution. We should mention that there have appeared many attempts to improve the 
rigid body approximation of [6,7]. They include the papers [9-29]. We wish to emphasize 
that our method based on the EOM principle gives a systematic and self-consistent way of 
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analyzing the spinning motion of a field theory soliton without any ad hoc physical pictures. 

Finally, we remark that the EOM principle has already been used for solitons in gauge 
theories even in the "rigid body approximation" before introducing the relativistic corrections. 
Namely, the dependence of the gauge field on the collective coordinates must be determined 
from the EOM of Aq, the Gauss law, which is non-trivial even in the rigid body approximation. 
The examples are the magnetic monopoles in gauge-Higgs systems (see [30] for a review) and 
the baryon solution [31,32] in the Sakai-Sugimoto model [33,34]. 

The organization of the rest of this paper is as follows. In Sec. 2, we summarize the rigid 
body approximation to the spinning Skyrmion and show how our EOM principle is violated. 
In Sec. 3, we introduce the improved spinning Skyrmion field with the leading relativistic 
correction to obtain, on the basis of the EOM principle, the differential equations and the 
boundary conditions of the functions (A(r), B(r),C(r)) specifying the correction. We also 
consider the deformation of the spinning Skyrmion caused by the relativistic correction. The 
Lagrangian of the spinning Skyrmion with the leading relativistic correction is given in Sec. 
4. In Sec. 5, we determine the values of the parameters (/„-, e) of the Skyrme theory from 
the masses of the nucleon, A and the pion. The profiles of the functions and the integrands 
related with the corrections are given. The deformations of the baryons are also pictorially 
shown there. Then, we present the analysis of the static properties of nucleons. The final 
section (Sec. 6) is devoted to a summary and future problems. In the Appendices (A-E), we 
present technical details and complicated equations used in the text. 

The present paper is a detailed version of [35], where the EOM principle was first proposed 
and applied to the quantization of the spinning motion of the Skyrmion. 

2 Skyrmion in the rigid body approximation and its 
limitation 

In this section, we summarize the Skyrmion and the quantization of its spinning degrees of 
freedom in the rigid body approximation [6,7]. We also explain how the EOM principle is 
violated in this approximation. This is the starting point of our relativistic extension. Another 
purpose of this section is to fix our notations and conventions. 

We consider the 577(2) Skyrme model [5] described by the chiral Lagrangian with the 
Skyrme term: 

C = tl {"IN + 32^' Lv]2 + f ml{U " l2) } ' (2 ' 1} 
where U (x) is an SU (2)-valued scalar field, is defined by 

L„ = -iUdptf, (2.2) 
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and e is a dimensionless parameter. Our flat space-time metric is r]^ = diag (— 1, 1, 1, 1), 
though not written explicitly in (2.1). As we will see later in Sec. 4, the inclusion of non-zero 
pion mass m n is indispensable for our relativistic corrections to be finite. The EOM of the 
Lagrangian (2.1) reads 

d, (V ~-^[Lu, [V, £']])- iml (u-±trU)=0. (2.3) 
The Skyrme model has a topologically conserved current J B : 



r B = -^- 2 e^tr(L u L p L a ), (2.4) 



with e 0123 = 1. The conserved charge is identified with the baryon number N B '- 

N B = [ d 3 xJ° B . (2.5) 



Next, the diagonal SU(2) symmetry of the theory gives the conserved isospin current: 

= tr [(■/£ + J" R ) r a ] , (2.6) 
where the left and right chiral currents J^/u are defined by 



(2.7) 



I L M replaced with _R M ' 

with = —iWdfJJ. On the other hand, the conservation of the axial 577(2) current, 

^, a = tr[(j£-j£)T a ], (2.8) 

is broken due to the pion mass term; d ll Jj ia = —i(f%m%/8) tr(C/r ). 

The Skyrme model has static soliton solutions which are called Skyrmion and represent 
the baryons [5]. The simplest one takes the hedgehog form: 

U c i(x) = exp(ix ■ rF(r)) = cosF(r) + ix ■ TsinF(r), (2.9) 

with r = \x\ and x = x/\x\. The function F(r) is subject to the following differential 
equation: 



d 2 F 2dF sin2F 4 
+ — ; ^ + 



dr 2 r dr r 2 f 2 e 2 



2 sin 2 Fd 2 F sin2F f dF\ 2 sin 2 Fsin2F 
+ s— ■ 



r 2 c?r 2 r 2 \ dr J r 4 



m 2 sinF = 0, 
(2.10) 
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which is obtained by substituting (2.9) into the EOM (2.3), or equivalently, by minimizing 
the mass (energy) M d of the configuration (2.9): 

nfl f°° , 2 f/i-i/\2 2sin 2 F 4 sin 2 F f sin 2 Fl 2/ 

M cI = ~ 2 ^ rfrr |( F ) + + [ 2(F + "^H + 2m ' (1 " COSF) j ' 

(2.11) 

with the prime on F denoting a differentiation with respect to r. For the hedgehog (2.9) to be 
non-singular both at the origin and the infinity and for the mass (2.11) to be finite, we must 
have F(0) = nn and F(oo) = (mod 2ir) with the integer n being equal to the baryon number 
N B . Restricting ourselves to the solution with unit baryon number N B = 1, the behavior of 
F(r) near the origin is given by 

F(r) = 7T - nr + 0(r 3 ) (r ->■ 0), (2.12) 

while that near the infinity reads 

F(r) = - ( 1 + — 1 e~ m - r + 0(e~ 2m " r ) (r -)• oo), (2.13) 

with k and a being constants. In particular, the power part multiplying e~ m7Tr in (2.13) is 
exact. The hedgehog solution (2.9) specified by a single function F(r) represents a spherically 
symmetric extended object; the energy density depends only on r. 

The hedgehog solution (2.9) has two kinds of collective coordinates; the center-of-mass 
motion and the space (isospace) rotation. Note that, for the hedgehog (2.9), the space rotation 
is equivalent to the isospace rotation. Namely, for any orthogonal matrix O, we have 

UciO-'x) = WU cl (x)W-\ (2.14) 

where W is the SU{2) matrix corresponding to O; O a bTb = WTaW -1 . In the quantization of 
the spinning collective coordinate in the rigid body approximation [6,7], we assume that the 
spinning Skyrmion field U(x,t) is simply given by 

U(x,t) = U cl (R- 1 (t)x), (2.15) 

where the orthogonal matrix R(t) is the dynamical variable to be quantized. Since the classical 
solution U c i(x) is of spherical shape, so does the Skyrme field of (2.15). An important property 
of U(x,t) (2.15) is that the left and the right SO (3) transformations on the matrix R(t) 
correspond to the rotations in the real space and the isospace, respectively: 

U M\ R ^o-l lR o iso = m so U(0^x : t)W.^. (2.16) 

Therefore, we write the components of R(t) by Ri a (t), with % and a being the real and the 
isospace indices, respectively. 
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The dynamics of R(t) is governed by its Lagrangian obtained by inserting (2.15) into the 
Skyrme model Lagrangian density (2.1) and carrying out the space integration. It turns out 
to be the Lagrangian of a spherical rigid rotor: 

L ligidhody (R,R) = J d 3 xC\ u{xA=UdiR _ lx) = -M cl + iztt 2 , (2.17) 



where 

2vr/2 roc 



dr r 2 sin 2 F < 1 



3 jo ~" ~~ ~ r 1 fy 

is the moment of inertia, and 



\ F 'f, sin2jP 



(2.18) 



V a = ^eaUR^R^c, (2.19) 

is the angular velocity in isospace (which is equal to the angular velocity in the body-fixed 
frame). The angular velocity Ui in real space is given by 

Ui = —-j^ijk{RR l )jk = &>i\ R ^ R -i — —Ria^a (2.20) 

Note that we have 

Q 2 = u 2 = —TriR^R) 2 . (2.21) 

The Lagrangian (2.17) describes baryons with unit baryon number N B = 1 and zero center- 
of-mass momentum. They are specified by the quantum numbers of the spin J = Xuj and the 
isospin i" = 2X2, which are the Noether charges corresponding to the symmetry transformations 
of the Lagrangian (2.17): 

R(t) O^R{t)O iso . (2.22) 

Since / and J are related by Jj = ^l^^-i = —Riala in the present system, we have I 2 = J 2 , 
which is consistent with the real baryon spectrum. The decay constant f n and various static 
properties of nucleons including charge radii, magnetic moments and axial vector coupling 
have been calculated in [6] (m n = case) and [7] by using the masses of the nucleon, A and 
the pion as inputs. The results agree with the experimental values within about 30% for most 
of the quantities. 

However, the quantization of the spinning collective coordinate explained above is not a 
fully satisfactory and consistent one as we mentioned in the Introduction. Here, we recapitulate 
the reasons: 

• Firstly, the Lagrangian (2.17) is nothing but that of a spherical rigid body. However, the 
rotating soliton should intuitively deform from the spherical shape, and the Lagrangian 
of R should contain terms which automatically incorporate such deformation. 
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By identifying the eigenvalues of the rigid body Hamiltonian (1.8) corresponding to the 
Lagrangian (2.17) with the experimental masses of the nucleon and A, we saw that 
about 8% (30%) of the total mass of the nucleon (A) comes from the rotational energy, 
and that the spinning velocity at the radius of 1 fm is nearly (or over) the light velocity. 
Therefore, the validity of the rigid body (i.e., non-relativistic) approximation for baryons 
is a non-trivial problem. 

The Skyrme field (2.15) does not satisfy the field theory EOM (2.3) even if we use the 
EOM of R(t) obtained from the Lagrangian (2.17): 



dfl d , • 

— = — R~ 1 R = 0, 
dt dt 



which is equivalent to 



du> d 



RR~ l = 0. 



dt dt 

In fact, the EOM (2.3) is violated by terms of 0(fl 2 ): 



(2.23) 



(2.24) 



LHS of (2.3)| 
2 



U(x,t)=U cl (y) \ ^« f 



+ 



+ i + 



X 



{F'f - 
sin2F 



2 

F " _|_ 
\ r 

1 - cos 2F 



1 - cos 2F) ( F" + -F' - + 2 sin 2F {F'f J ^{R^Ryfiy ■ r) 



2r 



{{R- l Rfy)r a 



1 - cos 2F 
2^2 



[y x {{R- l Rfy)]r a 



where y and L^(y) are defined by 



y = y(x,t) = R 1 (t)x, 
L*(y) = -iU d {y)-^U cX {y)\ 



(2.25) 

(2.26) 
(2.27) 



and we have used the EOM (2.3) for U c \. (For the derivation of (2.25), see Appendix 
A, where we present the derivation (3.17) with relativistic corrections.) Therefore, the 
quantization of spinning motion starting with (2.15) is valid only for slowly rotating 
baryons with small angular velocity fl. 



In the subsequent sections, we study relativistic corrections to the rigid body approximation 
on the basis of the EOM principle stated in the Introduction. 
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3 Realizing the EOM principle to O(d^) 



As we saw in the previous section, the spinning Skyrmion field given by (2.15) does not satisfy 
our EOM principle for introducing the collective coordinates: The field theory EOM (2.3) is 
broken by terms of 0(fl 2 ) = 0(8?) after using the EOM of R(t) (see (2.25)). Let us therefore 
try to resolve this 0(d 2 ) breaking of field theory EOM to shift it to the next higher order of 
0(df). Although we are still assuming that the rotation is not so fast, our analysis in this 
paper must prove an important step toward the complete realization of our EOM principle. 



3.1 Improved spinning Skyrmion field 

Our proposal for the improved Skyrmion field with the collective coordinate R(t) of spinning 
motion is that all the independence is contained in the argument of U c \ as in (2.15) in the 
rigid body approximation: 

U(x,t) = U cl (y(x,t)), (3.1) 

where y(x,t) is given by 

y(x, t) = \l + A^RR^xf + B{r)r 2 Tr(ir x i?) 2 + C (r)r 2 (R' 1 R) 2 ~\ R^x (3.2) 

with r = \x\. Compared with (2.15) in the rigid body approximation, we have introduced 
terms of 0(f2 2 ) in the argument of U c \. The present y(x,t) is a spinning motion analogue of 
(1.4) for the center-of-mass motion. The three functions A(r), B(r) and C(r) in (3.2) should 
be determined to fulfill our EOM principle. However, it is quite non-obvious at this stage 
whether the improvement can be accomplished by (3.1) and (3.2), namely, whether we can 
consistently determine (A(r), B(r),C(r)) so that our EOM principle is satisfied to 0(f2 2 ). 

The form of y(x,t) (3.2) is the most general one containing at most two time-derivatives 
and satisfying the following requirements: 

• y(x,t) is odd under x — > —x, 

y(-x,t) = -y(x,t), (3.3) 
and is even under the time- inversion of R(t): 

y( x ' t )\ R (t)^R(-t ) = y( x '- t )- ( 3 - 4 ) 

• y(x,t) has the following property under the constant left and right SO(3) transforma- 
tions of R{t): 

vMWt^o^mo-n = Or s oV{O rea iX,t) (3.5) 
This implies that our improved U(x,t) of (3.1) keeps the property (2.16). 
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• y(x,t) does not contain R(t). This is necessary for the Lagrangian of R(t) to consist of 
R and R without R. 

As we will see later in Sec. 4, the EOM of R(t) corresponding to the improved U(x : t) (3.1) 
remains the same as (2.23). Under the EOM (2.23), our Skyrmion field (3.1) is indeed spinning 
both in space and isospace with angular velocities u) (2.20) and Q (2.19), respectively. This 
fact may precisely be stated as follows. First, y(x,t) of (3.2) has the following property: 



§; + (RR-i X )°)y( X ,t) 



= 0, 



(3.6) 



R-EOM 



where "|_r_eom" means "upon using the EOM (2.23) of R(t)" . Since U(x,t) (3.1) is a function 
of y only, it obeys the same equation as (3.6): 



| + («x*),A)170M) 



0. 



(3.7) 



_R-EOM 



This implies that U(x,t) is spinning in real space with angular velocity u> when the EOM 
(2.23) holds. Next, y(x,t) also obeys 

= -R-'Ry = Vxy. (3.8) 

°^ _R-EOM 

This together with the hedgehog property (2.14) means that U(x,t) is spinning in isospace 
with angular velocity Q. 

Finally in this subsection, we consider the geometrical meaning of y(x,t) (3.2), and, in 
particular, the spatial shape represented by our spinning Skyrmion (3.1). Below we assume 
that r 2 Q 2 is sufficiently small. Let u be the displacement vector, namely, the difference of y 
(3.2) and y = R^x (2.26): 

y = y + u. (3.9) 

Taking, in the y-space, the polar coordinate system (r= \y\ = \x\ ,0,tp) with z-axis in the f2 
direction, u has the following decomposition in terms of the unit vectors e r = y/r and e e (see 
Fig. 1): 



u = 



Y(r) + (A(r) -C(r))P 2 (cos 6)je r + C(r) sin 9 cos 6 e e j r 3 fi 2 , (3.10) 

where the function Y(r) is defined by 

Y(r) = -A(r) + 3B(r) + C(r), (3.11) 
and P2 is a Legendre polynomial: 

P 2 (z) = \ (3z 2 - 1) . (3.12) 
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Note that the functions (A, B, C) appear in u as another three independent combinations; 
Y, A — C and C, which also appear in other places below. The decomposition (3.10) of u is 
restated as the following relation between the polar coordinate (\y \ ,9 y ,ip) of y and (r,0,tp) 
of y (the tp coordinate is common since there is no term in (3.10)): 



I y I 



1 - 



Y(r) + (A(r) - C(r))P 2 (cos9) 



r 2 ft 2 \ r, 



tan0 y = (l-C(r)r 2 ft 2 )tan0, 
where 9 y is the angle between y and fl. 




e 



(3.13) 
(3.14) 



Figure 1: Polar coordinate (r,9) for y, and the unit vectors (e r ,eg) in (3.10). 



The (classical) shape of our spinning Skyrmion (3.1) in a;-space for a given value of Q 2 
may have various definitions; for example, surfaces of constant energy density, constant baryon 
number density, and so onJ Here, noting that U c \{y) represents a spherically symmetric object 
in y-space, let us take the simplest one; surfaces of constant \y\. From (3.13), the surface of 
\y\ = a (= const.) is mapped to the a>space as 



5 = i + 

r 



(A(a) - C(a)) sin 2 9 - 2B{a) a 2 n 2 + 0{a 4 Q 4 ). 



(3.15) 



Here, we have r = \x\, and n — 9 is the angle between uj and x. Eq. (3.15) represents an 
ellipse on the (r,9) plane, and hence our spinning Skyrmion (3.1) has a shape of spheroid, as 
is intuitively expected. 



3.2 Differential equations and boundary conditions of (A,B,C) 

The three functions (A, B, C) in (3.2) must be determined in such a way that the field theory 
EOM (2.3) with the spinning Skyrmion field (3.1) substituted holds to 0(fl 2 ) upon using the 

^The isoscalar quadrupolc moment operator, which vanishes identically in the rigid body approximation, 
can be non-trivial due to present relativistic correction and is given by (E.38) in Appendix E. However, its 
nucleon expectation value is equal to zero. 
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EOM of the collective coordinate R(t). Namely, the contribution of the f2 2 terms in (3.2) to 
the field theory EOM must cancel the f2 2 terms in (2.25). 

After a tedious but straightforward calculation by using 

Lo(x,t) = ^^L^y), L t (x,t) = ^^^(y), (3.16) 

with L^(y) = —iU c \(y)(d/dy a )U c \(y)^ defined by (2.27), we get, instead of (2.25), the following 
result (we summarize the derivation in Appendix A): 

LHS of (2.3) | = (l? ~^[Lt [Ltm) M [(j t R " R ) V 
+ r 2 TriR-'Rfiy • r) x EQ X + {R^Ryfiy ■ r) x EQ 2 

+ r 2 [{R- 1 Rfy] ■ r x EQ 3 + r[y x [R^Rfy] - r x EQ 4 + 0(df), (3.17) 

where EQ n (n = 1,2,3,4) are linear in (A,B,C) and their first and second derivatives with 
respect to r with coefficients given in terms of F and its r-derivatives. The concrete expressions 
of EQ„ are very lengthy, and they are presented in Appendix B. As we have already mentioned 
in Sec. 3.1, and as we will see in the next section, the EOM of the collective coordinate 
R(t) derived from its Lagrangian, which is obtained by substituting (3.1) into the Skyrme 
Lagrangian density (2.1) and carrying out the space integration, remains unchanged from 
(2.23) in the rigid body approximation. This implies that the first term on the RHS of (3.17) 
vanishes upon using the EOM of R(t), and our EOM principle is satisfied to 0(fJ 2 ) if the 
three functions (A,B,C) satisfy four differential equations, 

EQ n = 0, (n = 1,2,3,4). (3.18) 

This is apparently overdetermined, but, fortunately, EQ 3 and EQ 4 are not independent. We 
have 

EQ 3 = 2cosF x EQ 34 , EQ 4 = -2 sin F x EQ 34 , (3.19) 

with EQ 34 given by (B.5). Therefore, (A, B, C) are determined by solving three inhomogeneous 
differential equations of second order, EQ X = EQ 2 = EQ 34 = 0. 

As we can guess from the fact that Y(r) (3.11) is the coefficient of the lowest mode 
1 = Po(cos#) in (3.13), we can extract the differential equation for a single function Y(r) by 
taking a suitable linear combination of EQ 12)3 . In fact, the special combination, 

EQ y = -3EQ 1 + EQ 2 -EQ 3 . (3.20) 

consists only of Y(r) and its derivatives as given in (B.6). Note also that B(r) and its 
derivatives are missing from EQ 2 (B.2) and EQ 34 (B.5). Therefore, for practical purposes, it 
is convenient to solve EQ y = for Y(r), and EQ 2 = EQ 34 = for (A(r), C(r)). 
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For solving the differential equations for (A,B,C), we have to specify the boundary con- 
ditions of the three functions at r = and r = oo. First, substituting the approximate 
expression of F(r) near the origin given by (2.12) into EQ 2 , EQ 34 and EQ y , we obtain the 
following differential equations for (A, C) and Y valid for r ~ 0: 

(1 + 8/i 2 )r 2 ^ - W 2 ^? + 8(1 + lf)r d 4- - 12A^ = 0, 
dr z dr z dr dr 

(1 + V)r 2 ^ - + 8(1 + 5/i 2 )r^ - 4(1 + 14/i 2 )A + (10 + 56/i 2 )C = 1 - 4/i 2 , 

with dimensionless ii defined by fx — n/ef^. From this we find that (A,C) and Y near the 
origin is generically given in terms of a particular solution and six independent modes: 

A(r)\ l-4fi 2 f 4/i 2 \ [5 + 28/i 2 \ f0\ 1 / -Q/j 2 \1 (5\ 1 
C(r)J " 10(1 + 8^) ^l + 4/i 2 y ' +Co ^2 + 28 / iV +C2 ^1 J r 2 + ° 5 [l + ty 2 ) r^° 7 [2)^ 

Y u = w^) + *h + *h> < r ~°>' (3 - 22) 

where c n (n = 0,2,5,7) and (n = 2,5) are constants. Each term on the RHSs of (3.22), 
for example, (|) (1/r 7 ), means a series in r 2 starting with this term. 

On the other hand, plugging the r — > 00 behavior of F(r) given by (2.13) into EQ 2 , EQ 34 
and EQ r , we obtain the differential equations for (A, C) and Y approximated near the infinity. 
They are lengthy and are given by (B.7) in Appendix B by using the dimensionless variable 
s = m^r. The Skyrme term does not contribute to these differential equations. From (B.7), 
we find that (A, C) and Y near the infinity are given as the following sum of a particular 
solution and six independent modes: 

y (0 = A + #i + tfi, (r^oo), (3.23) 



where rf 2j3 , / 2 3 , d\ and / 3 y are constants, and each term on the RHSs shows the first term of 
the series in 1/s. Concrete expressions keeping terms up to 1/s 5 for the power parts are given 
in (B.8) in Appendix B. 

For solving the differential equations for (A, B, C) globally in the range < r < 00, we 
impose the following boundary conditions in terms of the coefficients appearing in (3.22) and 
(3.23): 

c 5 = cl = c 7 = 0, (3.24) 
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h = h = H = o. 



(3.25) 



Namely, we choose the "mildest" boundary conditions both at the origin and the infinity. 
As we will see later, the most important functional of (A,B,C) is J given by (4.2). It is 
the coefficient of the f2 4 term of the Lagrangian (4.1) of R(t) and hence appears in various 
conserved charge operators such as Hamiltonian, angular momentum and isospin. Owing to 
our boundary conditions, (3.24) and (3.25), the r-integration for J is convergent at both ends, 
r = and r = oo. This is the case also for other physical quantities which we will discuss 
later in Sec. 5.3. Moreover, the deformation vector u(x,t) (3.10) is non-divergent at r = 0. It 
is also finite in the limit r — > oo if we fix r 2 Q 2 = finite. These properties look natural for our 
improvement to make sense. 

Having fixed the boundary conditions for (A,B,C), we can numerically solve the differ- 
ential equations for them by the shooting method. Namely, setting the boundary condition 
(3.24) at r = 0, we tune the three parameters (co, C2, c%) in (3.22) so that the condition (3.25) 
at r = oo is satisfied. This task of solving the differential equations for (A, B, C) will be done 
in Sec. 5.2 after we obtain the Lagrangian and the Hamiltonian for R(t). 

Finally, we add that our EOM principle implies that the spinning Skyrmion field (3.1) 
with a constant fl is nothing but a classical spinning solution of the Skyrme field theory up 
to 0(fi 4 ). 



4 Lagrangian of R(t) 



The Lagrangian of R(t) in our formalism with relativistic correction is given, as in the rigid 
body approximation of (2.17), by substituting the improved Skyrmion field (3.1) into the 
Lagrangian density (2.1) and carrying out the space integration. The corrections we have 
added in (3.2) contain two time-derivatives, and hence we can consider terms with four time- 
derivatives in the Lagrangian. We find that it is given by 



L{R,R) = d z xC 



\U(x,t)=U cl (y) 



-m c1 + ^m 2 + 



(4.1) 



where the rest mass M c \ and the moment-of-inertia X are the same as in the rigid body 
approximation, (2.17), and the coefficient J of the newly added f2 4 term is given by 



J 



15 

x 



dr r sin F 



rZ' + 5Z-C + 



sin 2 F 



(rZ' + 3Z + 2C) - {F'f (rZ' + Z + C) 



, (4-2) 
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with Z(r) defined by* 



Z = -2A + 5B + 2C = -Y - -(A - C). 

3 3 



(4.3) 



The Lagrangian (4.1) should be regarded as the first three terms of the rotational motion coun- 



terpart of the relativistic Lagrangian of the center-of-mass coordinate X(t); — M c \y \ — V 2 = 
-M c i + (l/2)M cl V 2 + (l/8)M cl V i + . . .. The EOM of R(t) derived from the Lagrangian (4.1) 



the rigid body approximation, which we already used deriving the differential equations for 
(A, B, C) in Sec. 3.2. 

Our result (4.1), in particular, J of (4.2), can be obtained by a rather lengthy calculation 
with use of (3.16) for L^. The outline of the derivation of (4.2) is given in Appendix C. Here, 
we explain a number of points used in deriving (4.1) and (4.2). 

First, note that Lagrangian (4.1) consists only of the angular velocity S7; R(t) and hence tl 
are missing. It is due to this fact that the EOM of R(t) remains essentially the same as in the 
rigid body approximation. The reason why the higher time-derivative of R(t) does not appear 
in the Lagrangian is as follows. Since the Lagrangian of R(t) has an invariance under (2.22) 
due to the property (3.5) and hence (2.16), and since y(x,t) (3.2) does not contain R nor 
terms with a single R _1 R, the possible terms in the Lagrangian of R(t) containing R and at 
most four time-derivatives in total are Tr^R^R (d/dt)R~ l R) and Tr ((iT 1 R) 2 (d/ dt) R' 1 R) . 
The origin of these terms is the part quadratic in Lq in the Skyrme model Lagrangian density 
(2.1). However, neither of these two terms can actually exist: The former term with three 
time-derivatives cannot appear in the Lagrangian owing to another property (3.4) of y(x,t) 
(this term is in any case a time-derivative term and can be discarded even if it does exist). 
The latter term vanishes identically since R~ X R is an anti-symmetric matrix. 

Now, the Lagrangian of R(t) takes symbolically the form L = 1 + fl 2 + fl 4 . To consider 
the coefficients of the three terms, we divide the Lagrangian into the "kinetic part" T (the 
part with L ) and the "potential part" V (the part without Lq), and write L = T — V with 



*The function Z(r) is extracted from the coefficient of e r in (3.10) by carrying out the cos ^-integration 



reads 





(4.5) 



with weight function sin 2 8: 
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V = Jd 3 x tr {|x 2 - ^L[ Li , Lj] 2 - £ml(U - 1 2 )} 



U(x,t)=U cl (y) 



(4.6) 



Then, T and V have the following expressions: 



and hence 



T = ±in 2 + ±j 1 n 4 + o(d?), (4.7) 

V = M cl - ^Am 2 - \j 2 ^ + 0(df), (4.8) 

L = -M cl + I (X + AX) ft 2 + i ( Ji + J" 2 ) ft 4 . (4.9) 

The origin of each term in (4.7) and (4.8) is explained as follows. Since y = R^x + u of 
(3.2) is the sum of the rigid body term R~ x x and the correction term u which are quadratic 
in R^R, the coefficients of the lowest order terms in (4.7) and (4.8), namely, X and M c \, are 
the same as those in the rigid body approximation. Next, the J7if2 4 and the AXfi 2 terms are 
the leading terms due to our relativistic correction u, and hence they are linear in (A, B, C). 
Finally, the ^T 2 Q A term in (4.8) is quadratic in u and therefore is quadratic in (A, B, C). One 
might wonder that, if we have added the 0(df) terms to y, they would also contribute to the 
f2 4 term in V. However, this is not the case as we will see below. 

By a simple scaling argument, we can get useful relationships among the coefficients in 
(4.9). Let us replace (A, B, C) in y of (3.2) by X(A, B, C) with A being a constant parameter, 
and denote this modified y by y\, namely, y\ = R^x + Xu. Then, using the fact that AX and 
Ji are linear in (A, B, C), while J2 is quadratic in them, the Lagrangian of R(t) corresponding 
to U(x,t) = U c \(yx(x,t)) reads, instead of (4.9), as follows: 

L\ = -M cl + A AX) Q 2 + - (A Ji + A 2 J 2 ) CI 4 . (4.10) 

At this point we should recall our EOM principle of introducing collective coordinates; U (x, t) = 
U c \(y\=i) satisfies the field theory EOM, namely, it is an extrema of the field theory action if 
R(t) is subject to its EOM, fl = 0. This implies, in particular, that L\ (4.10) is stationary at 
A = 1 for any constant fl: 



°=Tx L > 



= -Alfl 2 + -(J 1 + 2J 2 )fl\ (for v ft). (4.11) 

A=l 1 4 

Therefore, we get 

AX = 0, Ji + 2J 2 = 0. (4.12) 

By the same rescaling argument (by using another parameter), we can show that the addition 
of O(df) terms to u does not affect the coefficient of f2 4 in (4.9) as we mentioned above. § 



§The same kinds of relationships also hold in the case of the relativistic correction to the center-of-mass 
motion described by (1.4), (1.5) and (1.6): In the Lagrangian (1.5), C2 does not affect the coefficient of the 
X 2 term, the relation C2 + 2(— c|) = holds for the coefficient C2 — of X 4 (since we have C2 = 1/2), and C4 
docs not enter the coefficient of the X term. 
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Owing to the second relation of (4.12), the coefficient J of f2 4 in (4.9) can be expressed 
only in terms of J\. 

J = Ji + J 2 = \ji (4.13) 

This is very beneficial for obtaining J since J 2 , which is quadratic in (A, B, C), is harder to 
evaluate. In Appendix C, we present the derivation of J (4.2) by using (4.13). 

Some comments are in order concerning our results (4.1) and (4.2) and the boundary 
conditions (3.24) and (3.25) we have chosen for (A, B, C). First, from the boundary behaviors 
of F(r) given by (2.12) and (2.13), the r-integration of J (4.2) is approximated near r = and 
r = oo by f r ^ drr 6 [rZ' + 5Z - C + (2/t/e/ u ) 2 (2Z + C)\ and / r ~°°dr r 2 e - 2m - r (rZ' + 5Z -C), 
respectively. They are both convergent due to the boundary conditions c-j = at r = and 
(3.25) at r = oo. 

Our second comment is on the necessity of the non-zero pion mass. If the pion mass 
m n is zero, the exponential falloff (2.13) of F(r) near r = oo is changed to the power one 
F(r) ~ 1/r 2 . Correspondingly, choosing the mildest boundary condition at r = oo, we have 
the following leading behaviors near r = oo; (A,C) ~ (—3/4,-1/2) and Y ~ 1/4. This 
implies that the r-integration (4.2) for J is linearly divergent at r = oo in the case of m n = 0. 
This is the case also for relativistic corrections to other physical quantities which we discuss 
in Sec. 5.3. Therefore, the inclusion of non-zero pion mass is inevitable for our analysis in 
contrast to the case of the rigid body approximation [6]. (On the other hand, the qualitative 
behaviors of the various quantities near the origin r = are insensitive to m n .) 

Next, as described in Appendix C, the expression (4.2) for J has been obtained by changing 
the variables of integration in (4.1) from x to y, and carrying out the integration over the 
solid-angle of y. Therefore, r in (4.2) is the length of y; r = \y\. This way of switching to the 
y-integration is easier than to consider the original ^-integration. If we persist in carrying out 
the cc-integration, we obtain the expression of J which differs from (4.2) by a surface term: 



4tt f 2 

A t f =eq.(4.2)-— ^ 

I by cc-integration 1 v ' U 



, (4.14) 

r=0 



with r = \x\ on the RHS. From (2.12) and (2.13) for F(r), the surface terms of (4.14) at 
r = and r = oo are (up to constants) lim r ^ ^ 7 ^(^) and lim^oo r 3 e _2m " rr Z(r) , respectively. 
Fortunately, they both vanish safely owing to the boundary conditions (3.24) and (3.25) of 
(A,B,C). 

Our last comment is on AX in (4.9). Explicit calculation by using the y- integration and 
the EOM (2.10) of F(r) leads to the following expression of AX as a surface term: 



2sin 2 F 4 sin F ( sin 2 F 



r=0 
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= Ml L 2 (l- 4^t1 lim r 5 Y(r) + a 2 ml lim r 3 e" 2m ^F (r) 1 . (4.15) 
This also vanishes owing to the boundary conditions c( = = of Y(r). 

5 f n and the static properties of nucleons 

In this section, starting with the Lagrangian (4.1) of R(t) with the relativistic correction term 
(1/4)J7T2 4 , we first determine the parameters (f w ,e) of the Skyrme model from the masses of 
the lightest baryons, the nucleon and A, as well as m^. This task includes the determination 
of the functions F(r) and the newly introduced (A(r), B(r),C(r)). Then, we calculate various 
static properties of nucleons and compare them with the experimental values and also with 
the results of [7] without relativistic correction. 

5.1 Hamiltonian and spin/isospin operators 

Let us consider the standard canonical quantization of the system described by the Lagrangian 
(4.1) with dynamical variable R(t) G 5*0(3). For this, we take three independent variables 
£ A (t) (A = 1,2,3) (for example, the Euler angles) which parametrize the S0(3) matrix R(t); 
R = R(£ A (t)). Then, we have 

n a = I^ur-'r)^ = i A w Aa , n 2 = gAB {m A i B , (s.i) 

with w Aa and #ab(0 defined by 

/ ab 

The canonical momentum tta conjugate to £ A is 

n A = ^- = (l+jn 2 )g A B(Oi B , (5-3) 
and the Hamiltonian of the present system is given by 

H = n A i A -L = M cl + -Ifl 2 + ^jn 4 . (5.4) 

For identifying the eigenvalues of the Hamiltonian with the masses of the baryons, we have 
to relate Cl 2 with the isospin I and spin J, which are the Noether charges corresponding to 
the symmetry transformations (2.22). The isospin operator is given by 

i a = (i + jn 2 )n a , (5.5) 
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and the canonical commutation relation 7Tb(0] = i$B leads to 

{la, h] = i^abclc, [la, Rib] = itabcRic- (5.6) 

Similarly, the spin operator is 

Ji=(X + JQ 2 ) u t = -R ia I a , (5.7) 

and satisfies 

[Ji, Jj] = i€ij k I k , [Jj, Rja] = ieij k R ka , [I a , Jj] = 0. (5.8) 
Then, the equation relating fl 2 with the spin and the isospin is 

I 2 = J*=(l + jtf) 2 tf. (5.9) 

The Hamiltonian and the (iso)spin charges above have been obtained from the Lagrangian 
(4.1) of R(t). One may, however, wonder whether these conserved charges in the quantum 
mechanical system of R agree with those from the Skyrme field theory, namely, whether 
H, I a and Ji above agree with those obtained by substituting U(x,t) = U c \(y) (3.1) into 
the corresponding Noether charges from the Lagrangian density (2.1) of the Skyrme model. 
This agreement indeed holds owing to our EOM principle (see Appendix D for a proof). 
Furthermore, by using the equality of the isospins I a obtained in two ways, we can get the 
expression (4.2) of J directly from the calculation in the field theory side without relying on 
the relation (4.13). This is explained in Appendix E (see the paragraph below (E.13)). 



5.2 Determining (f n ,e) and (A,B,C) from (Mn, Ma, m n ) 

Now, let us consider determining (/,,-, e) from the Hamiltonian (5.4) together with the relation 
(5.9) by taking the masses of the nucleon, A and the pion, 

M N = 939 MeV, M A = 1232 MeV, m w = 138 MeV, (5.10) 

as inputs. Concrete procedure is as follows. We rewrite the differential equations for F and 
(A, B, C) and the integrations for X (2.18) and J (4.2) in terms of the dimensionless variable 
p = ef n r. Then, the differential equations and the integrations for the dimensionless quantities 
(M c i,X, J) defined by 

M cl = ^M cl , X=e 3 /.X, J = e>fU, (5.11) 

Jtt 

are characterized by a single parameter /3 = va^jiej^. For a given value of (3, we first solve 
the differential equation for F and then those for (A,B,C) to obtain (M c i,X, J7"). Next, by 
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identifying the Hamiltonian (5.4) for the nucleon (J 2 = 3/4) and A (I 2 = 15/4) with the 
inputs M N and Ma, respectively, we get the corresponding (/,r,e, m„-). Namely, we solve 



I 2 = (T + e^jh 2 ) 2 ^ 2 , (5.13) 
with dimensionless fl 2 defined by 

^ = ^ 2 , (5-14) 

for both the nucleon and A to obtain (/„., e) and m n = ef n (3. We tune (3 to reproduce the 
experimental value of m n . Our result obtained this way is 

/ 7r = 125MeV, e = 5.64, (5.15) 

which is realized at /3 — 0.196. Compared with the experimental value /i cxp ^ = 186 MeV, our 
f n is fairly improved from that of [7], fi AN ^ = 108 MeV, without relativistic correction. The 
dimensionless quantities corresponding to the result (5.15) are given by 

(M ch %J) = (37.9,70.2,279), (5.16) 

and the angular velocities of the nucleon and A are as follows: 

\n N \= 206 MeV (Q 2 N = 8.44 x 10~ 5 ), \tt A \ = 330 MeV (f2 A = 2.17 x 10~ 4 ). (5.17) 

Note that the angular velocity of the nucleon (A) has become larger (smaller) than that in 
the rigid body approximation (see eq. (1.9)). 

In the rest of this subsection, we present graphically various numerical results correspond- 
ing to (/ w ,e) of (5.15). 

Profiles of (A, B, C) 

First, in Figs. 2 and 3, the profiles of (A, B, C) as functions of ip — tan -1 ef n r (0 < ip < 
7r/2) are given. Fig. 2 shows the profiles of (A,B,C) themselves, and Fig. 3 shows those of 
fl 2 r 2 (C — A,Y,C), which are functions appearing in (3.10) and are directly relevant to the 
deformation of baryons. As Q 2 we have taken the nucleon angular velocity Q 2 N = (206 MeV) 2 . 



Profiles of the integrands of (M c i,X, J) 

Next, Figs. 4 and 5 show the profiles of the integrands of the dimensionless quantities 
(M c i,X, constituting the Hamiltonian. They are given as p-integrations; (2.11), (2.18) and 
(4.2) with the replacements r — > p, d/dr — > d/dp, — > (3 and (fn,e) — > (1, 1). In Fig. 4, 
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Figure 2: Profiles of A (red solid curve), B (black 
broken one) and C (blue dash-dotted one) at /3 = 
0.196. The horizontal coordinate is «/> = tan -1 p. As 
r — » — > 0), while A tends to a constant, B and 
C diverge as 0(l/r 2 ). As r -> oo -> 7r/2), all of 
them tend to zero. 



Figure 3: Profiles of (i) ft 2 N r 2 {C - A) (red solid 
curve), (ii) Q, 2 N r 2 Y (black broken one), and (iii) 
Q, 2 N r 2 C (blue dash-dotted one) at fi = 0.196 as 
functions of rp = tan -1 p. As r — > 0, the three 
functions all approach constants. As r — » oo, while 
J7^r 2 (C — A) and fi^r 2 !" go to the same constant, 
n 2 N /(2ml) = 1.11 (c.f., (B.8)), ft 2 N r 2 C diverges lin- 
early in r. 



we show the profiles of their integrands as the -^-integrations. Namely, giving O = M c \, X, 
(l/4)i7 as O = f^ 2 dip fo(ip)i Fig. 4 shows the profiles of the three foW (f° r the sake of 
visualization, we take (1/4) J7" instead of itself). On the other hand, expressing O differently 
as O = An f^ 2 dp p 2 go(?P), Fig. 5 shows the profiles of go(4>), which are the angle averaged 
densities in real space. We have foW — 47r(sin 2 if)/ cos A if)) go (4>)- 

Deformation of baryons 

As we explained in Sec. 3.1, the spinning Skyrmion field (3.1) describes a spheroid (3.15). 
In Fig. 6, we show the spheroids (ellipses) which represent the nucleon and A and correspond 
to the sphere of radius a = 1 fm without spinning motion. Namely, they are the spheroids 
(3.15) with fl 2 given by Q 2 N and Q\ of (5.17), respectively, and with a = lfm. In Fig. 7, we 
show the deformation vectors —u (3.10) of the nucleon representing the shift of each point on 
the original sphere due to the spinning motion at radii 0.5 fm and lfm. Note that y is on the 
sphere and y = y — u is on the spheroid. In both the figures, the angular velocity points to 
the vertical direction. 



5.3 Static properties of nucleons 

Having determined the parameters e) of the Skyrme model from the experimental values of 
(M N , Ma, m n ), let us calculate the various static properties of nucleons which were analyzed 
in the rigid body approximation in Refs. [6,7]. Concretely, we consider the charge radii, 
magnetic moments, magnetic charge radii, and axial vector coupling of nucleons. We are of 
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¥ 

Figure 4: Profiles of the integrands of tp- 
intcgrations for (i) M c \ (black dash-dotted curve), 
(ii) I (blue broken one), and (iii) (1/4)^7" (red solid 
one), at j3 = 0.196. Note that the profile for J is 
multiplied by (1/4). The horizontal coordinate is 
ip — tan^ 1 p. 
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Figure 5: Profiles of the angle averaged densities 
in real space of (i) M c \ (black dash-dotted curve), 
(ii) I (blue broken one), and (iii) (1/4) J (red solid 
one), at f3 = 0.196. The horizontal coordinate is 
if> = tan -1 p. 
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Figure 6: The spheroids (3.15) of the nucleon (red 
broken ellipse curve) and A (blue dash-dotted one) 
spinning around the vertical axis and corresponding 
to the original sphere (black solid circle) with radius 
a = 1 fin. The units on the horizontal and vertical 
axes are fm. 



Figure 7: Deformation vectors -u (3.10) of the 
nucleon at various angles 9. The red solid arrows 
and blue broken ones are — u on r — 0.5 fm and 
r = 1 fm, respectively. The units on the horizontal 
and vertical axes are fm. 
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course interested in how the relativistic correction modifies the numerical results from those 
of [7]. 

First, we need the analytic expressions of these static properties derived from the La- 
grangian density (2.1) of the Skyrme model and the spinning Skyrmion field (3.1). In Ap- 
pendix E, we present the definitions, outline of derivations and the final analytic expressions 
of the static properties. The analytic expressions are as follows: 

• Isoscalar mean square charge radius (r 2 ) /=0 : eq. (E.4). 

• Isovector mean square charge radius (r 2 ) I=1 : eq. (E.ll). 

• Isoscalar mean square magnetic radius (r 2 ) MI=0 : eq. (E.20). 

• Isovector mean square magnetic radius (r 2 ) M I=1 : eq. (E.28). This is in fact equal to the 
electric one (r 2 ) /=1 . 

• Isoscalar g-factor g I=0 : eq. (E.19). 

• Isovector g-factor g I=1 : eq. (E.26) together with (E.27). It has another and simpler 
expression (E.30). 

• Axial vector coupling constant qa'- eq. (E.34) together with (E.35) and (E.36). 

These quantities are all given as integrations over the dimensionless radial coordinate p = ef w r. 
The mean square charge radii are given in units of l/f 2 e 2 = (0.280 fm) 2 . In every quantity, 
the relativistic correction part is multiplied by e^fl% with being the dimensionless angular 
velocity of the nucleon (5.17). 

In Table 1, we summarize the numerical values of the static properties obtained by using 
(e, fifv) and the functions F and (A, B, C) determined in Sec. 5.2. Instead of the (/-factors, 
we present the nucleon magnetic moments ji p and \i n in units of Bohr magneton: 

Vp=\(9i=o + 9i=i), = ^ (gi=o - 9i=i) ■ (5.18) 

In the table, we also present the predictions of Ref. [7] without relativistic correction as well 
as the experimental values. 

As seen from the table, the difference between the prediction of this paper and that of 
Ref. [7] is in the range of 5% to 20%. (Note that each of our predictions is not given simply by 
adding the S7 2 correction to the result of Ref. [7] since the values of f n and e are also changed.) 
Although the fundamental parameter of the theory, /„., has been improved, the relativistic 
correction makes the value of theoretical prediction further away from the experimental one for 
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Prediction Prediction 



Experiment 



(this paper) (Ref. [7]) 



125 MeV 108 MeV 



186 MeV 




0.59 fm 0.68 fm 

1.17 fm 1.04 fm 

0.85 fm 0.95 fm 

1.17 fm 1.04 fm 
1.65 1.97 

-0.99 -1.24 
1.67 1.59 
0.58 0.65 



0.81 fm 
0.94 fm 
0.82 fm 
0.86 fm 
2.79 
-1.91 
1.46 
1.24 



9a 



Table 1: The static properties of nucleons. Prediction of this paper and that of Ref. [7] in the 
rigid body approximation both use the experimental values of (Mat, Ma, m n ) as inputs. We 
follow the notations of Ref. [6] . 

most of the static properties of nucleons (the only exception is (r 2 ) M/=0 ). However, this should 
not be regarded as a manifestation of problems of our basic EOM principle for introducing 
collective coordinates. There would be two possibilities for this unwelcome result. One is that, 
since the relativistic correction is rather large, we have to take into account the contributions 
from still higher orders in f2 2 for obtaining better results of the static properties of nucleons. 
Another possibility would be that we cannot expect the Skyrme theory, which is merely a 
low energy effective theory, to reproduce precisely the physics of the baryon sector even if 
the full relativistic treatment of rotational collective coordinate is carried out. Concerning 
the first possibility, we should note that the ratio of the contributions of the three terms of 
the Hamiltonian (5.4) to the baryon masses is 89 : 7 : 4 for the nucleon and 68 : 14 : 18 for 
A. This suggests us that the expansion in powers of the angular velocity Q is not a good 
approximation especially for A with a larger (iso)spin. 



In this paper, we proposed the EOM principle for constructing the relativistic system of 
collective coordinates of a field theory soliton, and applied it to the "next to the rigid body 
approximation" analysis of the spinning Skyrmion. The relativistic correction was successfully 
incorporated in the argument of the static Skyrmion solution in terms of the three functions 
(A(r) , B (r) , C (r)) obeying a consistent set of differential equations demanded by the EOM 
principle. We computed the decay constant and various static properties of nucleons with 
leading relativistic corrections by taking the masses of the nucleon, A and the pion as inputs. 



6 Summary and outlook 
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The relativistic corrections to the results in the rigid body approximation are found to be 
in the range of 5% - 20%. Though the value of the decay constant has become closer to 
the experimental one due to the correction, the results are not good for most of the static 
properties of nucleons. We also studied how the baryons deform from the spherical shape to 
the spheroidal one due to the spinning motion. 

Since we have introduced the relativistic corrections in a proper manner based on the 
EOM principle, our unwelcome result on the static properties of nucleons is an inevitable 
nature of the baryon sector of the Skyrme model: It should not be regarded as a problem 
of our basic EOM principle. If the simple Skyrme model can reproduce the real baryon 
physics more precisely, a remedy for the unwelcome results should be looked for in our use 
of expansion in powers of the angular velocity. Though we took into account only the first 
relativistic correction of 0(f2 4 ) in this paper, it would be necessary to consider higher order 
corrections, or to devise an approximation valid for an extreme large angular velocity, at least 
for analyzing A. For this, we have to consider a systematic generalization of the expressions 
(3.13) and (3.14) in terms of the Legendre polynomials. The analysis [8] of the spinning baby 
Skyrmion on the two-dimensional plane may be helpful for considering the extreme case in 
the 3 + 1 dimensional model. With a better (or the full) realization of the EOM principle, we 
can give more reliable predictions for the Skyrmion. It is also an important subject to apply 
our method to the collective coordinate quantization of other interesting physical systems. 
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A Derivation of eq. (3.17) 



In this appendix, we sketch the derivation of (3.17), namely, the LHS of the Skyrme field theory 
EOM (2.3) with the improved spinning Skyrmion field (3.1) substituted. In this derivation, 
the most important formula is (3.16) for L^(x,t) with L c ^{y) (a = 1,2,3) given by 

L*(y) = -iU cl {y)^-U cl {y)^ 

= -^[sin2F5 a6 + (2rF / -sin2F)y a y 6 -(l-cos2F)e a6c y c }r 6 , (A.l) 

where y a is a unit vector y a = y a / \y\. In this appendix, r denotes the length of y, and the 
argument r of the function F and its derivatives is omitted: 

r=\y\, F — F(r), F' = (A.2) 

As for the argument r of (A, B,C), the difference between r = \x\ and r — \y\ does not matter 
since (A,B,C) are already multiplied by the 0(d 2 ) quantities in (3.2). 

Necessary formulas 

Let us present the formulas of the differentiations of y(x,t) (3.2) with respect to x and t. 
First, the time- derivatives are given by 

dya< f ,1:) = -(R~ 1 Ry) a + 2 \A{R- l Ry) b {Wy) b + r 2 B Tr(WiT y a 
at L J 

+ r 2 C({ W, R~ 1 R}y) a + 0(d*), (A.3) 

= -(Wy) a + {{R- l Rfy) a + o(df), (A.4) 

where W is defined by 

w-ijr^A (A.5) 

Note that, in (A.3) and (A.4), we have presented their expressions in terms of y, which are 
obtained by rewriting the original ones in terms of y by using (3.2). The counting of the order 
of time-derivatives d t is of course is defined by regarding x (and not y) as a zero-th order 
quantity. Though we need only the first term on the RHS of (A.3) in this appendix, we have 
presented the expression valid to 0(df). This will be used in deriving the coefficient J of 
(4.2) in Appendix C. 

Next, for the a?- derivatives, we have 

djt ^ L = Ri b (*a b + T ab ) , (A.6) 
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d 2 y a (x,t) 



dxidxj 



(A.7) 



where T ab and T£ c are given by 
T ab = A l 

+ TriR^R) 2 



{R^Ryfda, - 2y a ((R- 1 R) 2 y) b \ + - — {Br x Ry) 2 y a y h 



Id A 



r 2 B 5 ab + [r 



— + 25 J y a y b 



+ (2C + r^j {{R-'R) 2 y) a y b + r 2 C {(R' 1 R) 2 ) ab + 0(df), (A.8) 



and 

-pa -pa u 1 ab r, a 

L bc — 1 cb= ~ — ~ Z/1 



dy c 



5 ab ((R- 1 R) 2 y) c + 5 ac ((R- 1 R) 2 y) b + ^((ir 1 ^) 2 ) 



be 



+ 



+ 



IdA 

r dr 
d 2 A 
dr 2 



(R-'Ry) 2 ^^ + 5 ac y b + y a (5 bc - y b y c )) -2y a (((i?- 1 i?) 2 2 /) 6 y c + ((R~ 1 R) 2 y) < 



Vb 



(R 1 Ry) 2 y a yby c 



+ 



(dB\ /3 dB d 2 B\ 

2B + r—\ Ti{R- l Rf{5 ab y c + 5 ac y b + y a 5 bc ) + ( -— + — J ^(R' 1 R) 2 y a y b y c 



{{R-'m^^m^b+iiR-'Rfy)^ 



In particular, we have 



(A.9) 



^ab + ^ba — 2 



AiR^Ryf + r 2 BTr(R- 1 R) 2 



Tab 



+ 2 



l^R-iRyf + (r^- + 2b) TriR-'R) 2 
r dr \ dr J 



y a yb 



- ( -2A + 2C + r^f- 
dr 



{(R^R) 2 y) a y b +{(R^R) 2 y) b y a \ + 2r 2 C{(R~ 1 R) 2 ) ab , 

(A.10) 



and 



9 d 2 B dB t \ , , 

2 ^ + 8r— -2A + KLBJ Tr^" 1 /?) 2 ^ 



+ ( r ^ + 8r f - 4A + 10C ) (^U- 
We also need formulas concerning L c ^(y). First, we have 
dLf(y) 



(A.ll) 



1 



-r 2 F" + (2 + cos 2F) rF' - - sin 2F 

2 



y a yby c 
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+ Q sin 2F - cos 2F rF'^j y a 5 bc + Q sin 2F - rF'^j (y b 5 ac + y c 5 ab ) 

+ (1 - cos 2F - sm2FrF')y a e bdc y d + X - (1 - cos2F) e a5c L. (A.12) 



Next, the double commutator of L^(y) is given by 
T abc (y) = -T acb (y) = [L*(y), [L?(y), L?(y)}} 



2(F'Y - 



1 - cos 2F 



y a |sin2F (y b 5 cd - y c 5 bd ) + (1 - cos2F) (y b e ced - y c e bed ) y e ^r d 



1 — cos2F f 

— < sin2F {5 ab 5 cd - 5 ac 5 bd ) + (2rF' - sin2F) (5 ab y c - 5 ac y b ) y d 

+ (1 - cos2F) {5 ab e ced - <5 ac e 6e d) & \r d . 



(A.13) 



Upon contraction, we have 

Thrift) : 



2 (F/)2+ l-cos2F 



sin 2F 5 ab + (1 - cos 2F) e acf ,y c I r b 



sin2F 



{F'f - 4 



1 - cos 2F sin 2F(1 - cos 2F) 
7, F H 5 



VaVbU. 



(A.14) 



Finally, the orthogonality, 

{R- 1 Ry) a y a = 0, (A.15) 

which is due to the anti-symmetric nature of the matrix R~ 1 R, often simplifies the calculations 
in many places. 

Evaluation of each term in (2.3) 

Let us obtain the expression of each term in (2.3) in terms of L a l (y) and its derivative. We 
keep only terms containing at most two time-derivatives. First, for d^L^ = —d L + dih^ we 
have 



d L (x,t) = ^L a \y) + ^^(d a L?)(y) 



dt 2 aKa/ dt dt 
-(Wy) a + ((R-'Rfy)^ L a \y) + (R~ 1 Ry) a (R~ 1 Ry) b (d a Lf)(y), (A.16) 

d 2 y a 

dxf aya/ ' dxidxi 



dM^t) = ^L a \y) + ^^(d a L?)(y) 



= (d a L*)(y) + Vt b L a \y) + (T ab + T ba ) (d a L?)(y), 
with (d a Lf)(y) = (d/dy a )Lf(y). 



(A.17) 
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Next, let us consider the part of EOM coming from the Skyrme term: 

d, [L v , [L", L u ]\ = -do [L u [L , Li]] - di [L , [L t , L }] + d t [L v [L t , L 3 \] . 



(A.18) 



For the first two terms on the RHS, we can substitute the lowest order expressions Lq 
-(R- 1 Ry) a L^(y) and L t = R ia L*(y) to get 

d [L h [L , Li]] = [-(Wy) a + {(R^Rfy)^ T bab (y) + {R- 1 Ry) a {R- 1 Ry) c -^T bab {y), 



di[L , [Li,Lo]] = — {{R- 1 Ry) a {R- 1 Ry) c T abc {y)) , 

with W defined by (A.5). The last term (A.18) is given, using (A. 6) and (A.7), by 

d ( dy a dy b dy c 



(A.19) 
(A.20) 



d 



T ba b + Fc C T bab + T c ab T abc + (T ac + T ca ) 



dy a 

where we have used T bc T abc = 0. 
Total of the LHS of (2.3) 

From the above results, (A.16)-(A.21), we finally obtain 

LHS of (2.3)| = -^-im, f U*(v) "trl/ d (») 



f d d 
I - — T bab + — — T abc 



\dy c 



dy b 



(A.21) 



dy a 



+ 



(Wy) a - {{R- 1 Rfy) a + T a bb ) C*(y) + {^—{R~ 1 Ry) a {R~ 1 Ry) b + T ab + T 6a ) 

dT abc (y) d 
dy b dy, 



with defined by 



T c ab T a Uy) + (T ac + T ca ) 



{R- l Ry) a {R- l Ry) c T abc {y) 



rcl ±_ 

a fy 



Tj 



bab- 



The first two terms containing no time- derivatives on the RHS of (A. 22) vanish: 

dy a 



im w [ U d (y) - -tr U d (y) ) = 0. 



dy b 
(A.22) 
(A.23) 

(A.24) 



This is nothing but the EOM of U c \. The next term, (yVy) a Cf, gives the first term on the RHS 
of (3.17), and the remaining terms give the EOM-breaking part containing EQ n (n = 1, 2, 3, 4). 
In the latter calculation, the orthogonality (A. 15) greatly simplifies our lengthy task. 

Finally, by dropping the terms containing (A,B,C) in (A.22), namely, by removing T ab 
and T bc , we obtain (2.25) before introducing the improvement. 
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B Expressions of EQ n=1 2 ,3,4 



In this appendix, we present concrete (and lengthy) expressions of the four quantities EQ n (n = 
1,2,3,4) appearing in (3.17). We also present the approximate expressions of the differential 
equations for (A, C) and Y near the infinity and their solutions. The primes on F denote 
differentiations with respect to r. First, EQ n are given as follows: 

EQ X = -F'A - F'4? — 2 ( F" + -F') ^---( 3F" + -F' - 4 sin2F^| B 



r 2 dr 2 \ r ) dr r \ r r 2 

^ ( 2F ' ~ \ ™ 2F ) C + i| ^ - ^s2F)F'A - 1(1 - cos2F)F'^ 



8 



f 2 e 2 ^ r 4 

1 - cos 2F) ^F" + ^F'^j + sin 2F (F' 
1 - cos 2F) ^4F" + ^-F' - 1 sin 2F^ + 4 sin 2F (F' f 



dB 
dr 



B 



_i (1 _ cos2 ^_± 

r^ dr r* 



(1 - cos 2F) IF" + -F' — \ sin 2F ) + sin 2F {F'f 
\ r r 2 J 



(B.l) 



EQ 2 



dr A V r dr r 



-3F" + -(cos 2F - 3)F' + — sin 2F 



+ 



3F" + -(1 -cos2F)F' 
r 



2F" + -(7 - cos 2F)F' - — sin 2F 



C 



r* 
dC_ 
dr 



A 



+ 



1_ | l s [(i _ cos 2F) (f» + V - A sin 2f) + 2 sin 2F (F' 

- 4(1 -cos2F)F'44 - 4 [2(1 - cos 2F) ( F" + -F'] + sin2F(F 
r l dr 1 r z \_ \ r J 

-1 (l-cos2F) ^8F" + -^(1 -2cos2F)F'- -|sin2F^ + 5sin2F(F') 2 
+ 1(1 - cos 2F) ^4F' - isin2F^ 



dA 
dr 



d 2 C 
dr 2 



+ 



4 

+ -5 



(l-cos2F) (4F" + i(7-2cos2F)F / - 4sin2F ) + 2sin2F(F') 2 
\ r r l ) 



(l-cos2F) 5F"--(l + cos2F)F' + — sin2F + 4sin2F(F 



,V\2 



EQ 3 = 2 cos F x EQ 34 , 
EQ 4 = -2sinF x EQ 34 , 



c), (B.2) 

(B.3) 
(B.4) 
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where EQ 34 in (B.3) and (B.4) is given by 



12 1 d 2 C 1 / 3 \ dC 

EQ 34 = — sinF + — cosFF'A sinF— — - - cosFF' + - sinF — 

2 r -i r z 2r dr 2 r \ r I dr 



1 / i \ i |2 

- — (4 cos F F' + - sin F j C + — sin F I ~", 

+ -JF'f^ + -F' Uf 1 + - sm2F'] A 
r 2 dr r A \ r J 



(F'Y - -(1 - cos2F) 



(F'y + - sin 2F F' + — (1 - cos 2F) 



dC_ 

dr 



4 1 
(F'Y + - sin 2F F' - — (1 - cos 2F) 



C . 



(B.5) 



Next, the combination EQ y (3.20) is given by 



EQ y = -3EQ X + EQ 2 - EQ 3 
\ , 8 sin 2 F^ 



/ 2 e 2 



+ \ IF" + -F' + — 

r fie 2 



dr 2 
sin 2 F 



2sin2F 16 
; — + 



F" + -F') + 

r 

2 



3 

r 

sin 2 F 



2 sin 2 F ( ,. 2 . 
— F + -F 



sin 2F 



/\2 



-(F') 



/ 2 e 2 



dr 
sin2F 



4F" + -F' 

r 



sin2F\ sin2F / 



2sin2F 



(F'Y 



Y 
(B.6) 



By using (2.13) for F(r), the differential equations EQ 2 = EQ 34 = for (A,C) and 
EQ r = for Y are approximated near the infinity (r — > oo) as follows: 



2 2\d 2 A n / 1 2 2\dA /3 7 12 12'. , 
1 + - + -T —T + 2 -1 + - + - + - — -2 - + - + — + — )A 



s s 2 ) ds 2 



s s" 



ds 



s s 2 s 3 s 4 



3 3 \ d 2 C ( 3 3 \ dC (I 3 6 6 , 
1 + - + -] — + 2(1----] — + 6 - + - + - + - C = 0, 



s s z 



ds 2 



s 2 s 3 I ds 



s s 2 s 3 s 4 



,1 2 2 \ A ( l\d 2 C ( 1 1 \ dC (A 7 7 . 
4 - + ^ + ^U- l + - ^r + 21 ~ ^ + 2 - + ^ + ^ C 



s s 2 s 3 J V s) ds 2 V s s 2 1 ds \s s 2 s 3 

1 1 
+ — + — = 0, 



s 2 s 3 
2 2 \ d 2 Y o 

s s 2 ) ds 2 \ s s° 



2 \ dY 



3 2 



1 ' s 3 ) ds \ s^ s 2 ^ s 3 ^ s A ) ^ s 3 ^ s A 



1 1 



(B.7) 



32 



where s is the dimensionless variable s = m n r. The general solution to (B.7) is given symbol- 
ically by (3.23) as the sum of a particular solution and six independent modes. More concrete 
expression keeping terms up to 1/ s 5 for the power parts is follows: 

(A(r)\__±( 1 + 3/s 3 1 - 3 / s2 Urfl|" 1 + 2 ^ 

\C(r)J ~ 2s \1 - l/s + 3/s 2 + 3/s 3 ) + d2 s 2 [i - 2 /s - A/s 2 ) + ^ s 3 ^2 + A/s ) 

+ ! f L( i-i7/(3 S 2 ) V f If3-14MK 

+ \ /2 s 2 _ 4 /(3s) - 2/(3s 2 ) y ) +/3 S 3 ^2-8/ Sy l J ' 



C Derivation of eq. (4.2) for J" 



In this appendix, we outline the derivation of eq. (4.2) for J , the coefficient of the f2 4 term 
in the Lagrangian of R(t). Here, as in Appendix A, r denotes the length of y; r — \y\. As 
explained in Sec. 4, calculation of J is reduced via (4.13) to that of J7i, the part of J linear in 
(A,B,C). It is convenient to use the expression of the Lagrangian L as the integration over 
y instead of x: 

L= [d 3 yj£ , (C.l) 



(C.2) 



U( X ,t)=U cl (y) 

where J is the Jacobian and is given via (A. 6) in terms of T ab : 



J 



Explicitly, we have 



det ( 


-1 


(a,i) \dXiJ 





det(5 ab + T ab ) = l-TrT + 0(«9 4 ). 

(a,6) 



Tr T = T aa = [hA -2C + r^-(A-C^J (R^Ryf+r 2 [hB + C + r^j Tt(R^R) 2 (C.3) 



Let us identify the part of JC in (C.l) which is quartic in f2 and linear in (A,B,C) (we 
call such part "QL-part" hereafter). We employ (3.16) for together with (A. 3) and (A. 6) 
for dy a /dt and dy a /dxi, respectively. First, for the trL 2 part of C (2.1), we have 



Jtr L\ 



JtrZ4 



QL-part 



- Tr T (R~ 1 Ry) a (R~ 1 Ry) b tr{L C a{y)Lf(y)), 



QL-part 



= 0. 



(C.4) 
(C.5) 



In obtaining (C.4), we have used the fact that the terms in (A. 3) containing W (A. 5) do not 
contribute to the Lagrangian (C.l), which we showed in Sec. 4. This can of course be confirmed 
by explicit calculation. For (C.5), we have used that the 0(<9 4 ) term of T ab) which is not given 
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in (A. 8), is quadratic (A,B,C) and hence cannot contribute to the QL-part. Similarly, for 
the Skyrme-term tr[L M , L u } 2 = — 2tr[L , Li] 2 + tr[Lj, Lj] 2 , we have 



J tr[L , Li 



J tr[_Z^, Lj 



QL-part 



QL-part 



- 8 cd TrT — (T cd + T dc )^(R~ 1 Ry) a (R~ 1 Ry) b tr( < [L a \Lf][Ll\L d l ])(y), 

(C.6) 

0. (C.7) 



Then, the QL-part of the Lagrangian (C.l) is obtained by using (C.4)-(C.7), (A. 8), (C.3) 



and 



i _ pns 2 F 

tr(L^) (y) = 2(F') 2 y a y b + (5 ab - y a y b ) , 



(C.8) 



tr ([L* Lf] [Lf, L c d 1 ]) (y) = -2 1 ™ 2F 5 ah (V) 2 ^ + 1 ™ 2F (5 cd - y c y d ) 

n (I -cos2F\ 2 p p , . 

+ 21 — 1 d ad d bc + (terms containing y a and/or y 6 ) , 

(C.9) 

and carrying out the solid-angle integration of y. Denoting by overline the solid-angle average, 

O = i- Jdtt y O, (CIO) 

the following formulas are of use: 

(C.ll) 
(C.12) 
In particular, we have 

(C.13) 



VaVb = y^afc, 

r 4 

VaVbVcVd = 77 (<5a6<5cd + S ac S bd + 5 ad 5 bc ) 

10 



(R-iRy) 2 = -jTriR-'R) 2 = jW, 



(R-iRy) 4 
where we have used 



r 
15 



TriR^R) 2 ^ 2 + 2 TiiR^Rf 



15 r iZ ' 



Tr(R- L RY = -2fT, 

(R^R) 3 = 1-Tt(R- 1 R) 2 R- 1 R. 



(C.14) 

(C.15) 
(C.16) 



By identifying the QL-part of (C.l) with (l^Jifi 4 , we obtain J = (1/2) Ji as given by (4.2). 
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D Equivalence of charges from L(R,R) and from the 
Skyrme field theory 



In this Appendix, we present a proof of the equivalence of the conserved charges (especially, the 
isospin I and the Hamiltonian H) obtained in two different ways; one is from the Lagrangian 
L(R,R) (4.1) of R(t), and the other is by the substitution of U(x,t) = U c \(y) (3.1) into the 
Noether charges from the Lagrangian density (2.1) of the Skyrme field theory. We show that 
the two conserved charges agree with each other up to the EOM of R(t) and that in the 
Skyrme field theory; the latter also vanishes upon using the EOM of R(t) due to our principle. 

Isospin charge I a 

First, let us consider the isospin charge I a . For this, recall that the Lagrangian of R is 
related to the Lagrangian density C(U(x,t)) of the Skyrme model by 



L(R,R) = I d 3 xC\ 



U(x,t)=U cl (y)- 



(D.l) 



The isospin charge I a on the i?-system side is obtained by considering the following time- 
dependent transformation 5\ acting on L: 

5 x R(t) = R(t)X(t), 5 x R-\t) = -X^R-^t), (D.2) 

where X a b{t) is an infinitesimal anti-symmetric matrix. We have 

8 x L(R,R) = X a (t)I a , (D.3) 

with X a {t) = {l/2)e abc X bc {t). 

For relating the isospin charge of (D.3) to the one on the field theory side, let us consider 
the response of C(U(x,t) = U c \(y)) under the transformation 5\ of (D.2). Since we have 
S^R^x) = —XR~ l x and 5\(R^ 1 R) = X + A], the transformation of y (3.2) consists 

of two parts: 



S\y = ~Xy + X ab Z ab (x, t), 
where Z ab is a three dimensional vector for each (a, b): 



(DA) 



[Z ab ) c = A(R~ 1 Ry) a y b — Br 2 (RT 1 R) ab 



5 ac (R l Ry) h + (R 1 R) ca y b \ -(a**b), 

(D.5) 

with y = R~ l x. The first part S^y = —Xy induces the standard isospin transformation on 
U(x,t) = U cl (y): 

(D.6) 



5f ] U cl (y)=i U cl (y),X a (t)^ 
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and therefore we have 



6®£(U{x,t) = U*{y))=\ a jfi > 



(D.7) 

U(x,t)=U cl (y) 



where Jy a is the isovector current (2.6) of the Skyrme field theory. (If A a in (D.6) depends on 
both t and x, we have 5 { ° ] C = d^\ a (x,t) J$ a instead of (D.7).) 

The second part of the transformation (D.4), S^y = \ a bZ ab (x, t), needs careful treatment. 
Let us divide the Lagrangian density (2.1) into two parts, C = C (L^) + (f 2 /8)m 2 tr (U — 1 2 ), 
and consider first the transformation of C Q {L^) consisting only of = —iU c \{y){d/dx ,x )Ul l {y). 
Using 

S^U cl (y) = \ ab (Z ab ) c -^-U cl (y), (D.8) 



we have 

with 77 defined by 
This gives 

with 



S?L^^ + i[L,(x,t),r)}, (D.9) 
77 = \ ab (Z ab ) c L?(y). (D.10) 

g = "f (»-±[L.,V,n). (D.12) 

In 5^ J (fix Co, we carry out the integration-by-parts for the fx — 1,2,3 parts of (D.ll). For 
the fx = part, notice that d/dt acting on (Z ab ) c Lf(y) in (D.10), which consists of y and 
R _1 R, can be expressed as 

— = -(RR^x)i— + (R-EOM terms), (D.13) 

where (i?-EOM terms) denote terms which vanish upon use of the EOM (2.23) of R (c.f., 
eq. (3.6)). Since (D.13) holds also against dC /dL and since we have (RR -1 )^ = 0, space- 
integration- by-parts for all fx in (D.ll) is effectively allowed to give 

5^ J d 3 xC = (A-term) - j d 3 x tr^ ^df" ) + (#-EOM terms). (D.14) 

The A-term is another contribution from the fx = term of (D.ll); d/dt acting on X ab in 77 
(D.10). However, this A-term in fact vanishes. This is because, after carrying out the solid- 
angle integration of y, the possible quantities multiplying anti-symmetric matrix X ab are all 
symmetric ones, 5 ab Ti{R~ l R) 2 and ((R^R) 2 )^, in our present approximation of keeping at 
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most four time- derivatives in the Lagrangian (see the argument given in Sec. 4 for excluding 
(d/dt)(R- l R) from the La grangian of R ). 

For the 5^ transformation of the remaining pion mass term in C (2.1), we have, using 
(D.8), 

5?tr(U- 1 2 ) = \ ab (Z ab ) c tr^^± = -itr (rjU cl (y)) . (D.15) 
From this and (D.14), we have the following for the whole C: 



6^ Jd 3 xC = -Jd 3 x tr 



7/ 



d dC fl 2 ( 1 

+ i—m^ U — - tr U 



dx» dL^ 8 ff V 2 



+ (i?-EOM terms). 

(D.16) 

Note that the quantity inside the square bracket multiplying 77 in the first term on the RHS is 
nothing but the EOM (2.3) of the Skyrme field theory. Due to our EOM principle, this term 
also vanishes upon use of the EOM of R. 

We have calculated the 5\ transformation of the both hand sides of (D.l); the transforma- 
tion of the LHS is (D.3), while that of the RHS is the sum of (the space-integration of) (D.7) 
and (D.16). An important point is that the whole of (D.16) vanishes when the EOM (2.23) 
of R holds. Therefore, comparing the coefficients of A a (t), we find that the isospin charge I a 
obtained from the Lagrangian (4.1) is equal to that in the Skyrme field theory up to the EOM 
of R: 

I a = fd 3 xJ^=°(x,t) + (R-EOM terms). (D.17) 

J ' U(x,t)=U cl (y) 

For the spin Jj, we consider the infinitesimal left transformation on R(t) instead of the 
right one (D.2). Then the equivalence of Jj between the R system and the Skyrme field theory 
can be shown quite similarly. 

Hamiltonian H 

Next, we show that the Hamiltonian (5.4) of the system of R is equal to the Hamiltonian 
of the Skyrme field theory with the substitution (3.1) up to the EOM of R. The proof is quite 
similar to the above one for the isospin, and we will describe only the points of the proof. 

By considering the time translation with infinitesimal time-dependent parameter e(t), 

5 £ R{t) = e(t)R(t), (D.18) 

the Hamiltonian (5.4) of the R system is obtained as 

5 £ L = eH + d t (eL). (D.19) 

To relate this H to the Hamiltonian of the Skyrme field theory, we consider the action of the 
same transformation (D.18) on the Lagrangian density on the RHS of (D.l). We have 

S e y(x, t) = e(t)d t y(x, t) + s(t)z(x, t), (D.20) 
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with 

z = 2 \A{R^Ry) 2 + Br 2 Ti^R) 2 + C^iR^R) 2 ] y. (D.21) 
The first part of the transformation, di^y = ed t y, gives simply 

5^ [ d 3 xC(U(x,t) = U cl (y)) = eH Skyrmc + d t (eL), (D.22) 

J U(x,t)=U cl (y) 

with i?skyrme being the Hamiltonian of the Skyrme field theory. For the second part of the 
transformation, 5^y = ez, the same equation as (D.9) holds by replacing rj (D.10) with 
ez a L^(y). Then, the same argument by using (D.13) applies also to leading again to 
(D.16). In this case, the absence of the e term is owing to that the possible quantities mul- 
tiplying e are Tr(i? _1 i?) 3 and Tr(i? _1 i?) 2 Tr^" 1 ^), which however all vanish identically. In 
this way, we obtain the equivalence of H and i?skyrme|c(a:,t)=i7 c i(i/)- 

The equivalence of the Hamiltonians obtained by two ways, which holds owing to the 
EOM principle, provides us with another and simple derivation of the relation (4.12). First, 
the Hamiltonian H corresponding to the Lagrangian (4.9) is given by 

H = M cl + i (X + AX) fl 2 + ^ ( Ji + J 2 ) fl 4 . (D.23) 

On the other hand, the Hamiltonian of the Skyrme field theory is obtained from the Lagrangian 
by inverting the sign of its "potential term" V (4.8): 



= T + V = Mci + I (X - AX) Q 2 + - ( Ji - J 2 ) n 4 . (D.24) 

U( X ,t)=U cl (y) 2 4 



-^Skyrme 

Comparing the above two, we reobtain (4.12). 



E Analytic expressions of the static properties of nucle- 
ons 

In this appendix, we outline the derivation of the analytic expressions of the various static 
properties of the nucleons. Throughout this appendix, r denotes the length of y, r = \y\, 
except (r 2 ) (= (x 2 )) denoting the charge radii of various kinds. 

Iso scalar charge radius 

The isoscalar mean square charge radius (r 2 ) I=0 is given in terms of the baryon number 
density J^(x,t) and a nucleon state \N) by 

(r 2 ) I=Q = (N\ Jd*xx 2 J B (x,t)\N}. (E.l) 
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From (2.4) and (3.16), we have 



J° B (x,t) 



24ir 2 



e ijk tr (L i (x,t)L j (x,t)L k (x,t)) 



(a,i) \ (JXi J 



J c B°(\y\), 



where J B 10 is the baryon number density of the static solution (2.9): 

1 l-cos2FrfF(r) 



(E.2) 



(E.3) 



47T 2 r 2 dr 

Owing to the presence of the Jacobian in (E.2), it is convenient to switch to the ^/-integration 
to evaluate (E.l) as follows: 



<r 2 ) /=0 = (iV| J<Pyx 2 J* (\y\)\N) 



= 4tt / drr A (N\ l — 2(A — C) (R~ 1 Ry) 2 — 2Br 2 Tr(i? _1 i?) 2 \N) J c B l0 (r) 



4ir 



where we have used the relation 



dpp*[l + -eW N p 2 Y) jr( P ) 



X 



1 - 2 (A - C) (R^Ry) 2 - 2Br 2 Tr^R) 2 



\y\), 



(E.4) 



(E.5) 



valid to 0(d 2 ), and (C.13) for the angle average over y. The final expression is given in terms 
of the dimensionless p = ef w r, Q 2 N (5.17) and the function Y(r) (3.11). 

Isovector charge radius 

To obtain the expression of the isovector charge radius, we first need that of the isospin 
charge density Jy a of (2.6) with J£, R defined by (2.7). Corresponding to (3.16) for L^, 
R^ = —iU^dfJJ for our spinning Skyrmion field (3.1) is given by 



Ro(x,t) 



dy a (x,t) d . . dy a (x,t) d 

- R a (y). tufa t) = — 5 R a (y)> 



with 



dt 



d 



(E.6) 



R c Xy) = -iU cl ( y y—U cl (y) 

= ^-{sin2F<5 ab + (2rF' - sin2F) y a y b + (1 - cos2F) e abc y c }r b 
= -L*(-y). 



(E.7) 



Note that U c \{y)^ = U c \(—y). Since we are considering the on-shell R(t) with W = 0, (A. 3) is 
reduced to dy a (x,t)/dt = —(R~ 1 Ry) a in the present calculation. This leads to 



L Q + Rq = i (1 - cos2F) e abc (R 1 Ry) a ybT c . 



(E.8) 
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As for the part in J v from the Skyrme term, we have 

[L it [L , L^] + (L M = - (5 ac + T ac + T ca ) (R- l Ry) b T abc + (y -> -y), (E.9) 

with T ac + T ca and T abc given by (A. 10) and (A. 13), respectively. From (E.8) and (E.9), we 
obtain 



16e 2 J v (x,t) = f 2 e 2 (L + R ) - ([£,, [L , L,]] + [i*, [i? ,^]]) 



= 2 



1 - cos 2F 1 



+ 



+ 



d 



1 - cos2F 



4 ( 3(A - C) + r— (A - C)J (F ) + 2(4 + C) 
I ( 35 + r^pj r 2 (F') 2 + (2B + C)(l - cos2F) 



1 - cos2F 



(R-'Ryf 



TriR^R) 2 }e abc (R~ l Ry) a y b r c . 

(E.10) 



Using this, we get the following expression for the isovector mean square charge radius of the 
nucleon: 



(N\fd*xx 2 J v>a=3 (x,t) \N) 



3^ fy Jo 

2 



(N\I 3 \N) 
dpp A sm 2 F{ l + 4(F') 2 + 4 



sin 2 F 



+ Z /tf N p 2 



sin 2 F 



P Z' + 7Z-C + 4{F'y {-pZ 1 + Z-C)+ 4 — (pZ' + 5Z + 2C) 



(E.11) 



with prime denoting d/dp and Z defined by (4.3). In deriving (E.ll), it is convenient to switch 
to the ^-integration by using d 3 x = (1 — TrT) d 3 y with TrT given by (C.3), and use (E.5) 
for x 2 and the following formulas for the angle averaging: 



e a bc (R- 1 Ry) b y c = -re abc (R l R) bc = -rfi a , 



e abc (R- l Ryf(R~ l Ry) b y c = --r 3 Ti(R- 1 R) 2 e abc {R~ l R) bc = -r 3 n 2 Q a 



(E.12) 



In rewriting Q3 in the numerator of (E.ll) in terms of ^3, we have used (5.13) for the nucleon, 
namely, 



(i + jn 2 N ) 



2 A" 1 



\n 



N 



(E.13) 



Here, we comment on another derivation of (4.2) for J . In Sec. 5, the isospin charge I a 
(5.5) was obtained from the Lagrangian (4.1) of R(t). As explained in Appendix D, the same 
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I a should also be obtained by integrating (E.10) over x. By the latter way, we directly obtain 
(4.2) for J (and also (2.18) for X) without using the relation (4.13). 

Isoscalar magnetic moment and magnetic charge radius 



Let us consider the isoscalar magnetic moment defined by 

M/=o = ^ Jd 3 xx x J B (x,t). 
The space component of the baryon number current (2.4) is given by 



(E.14) 



40M) = -^e^MLjLkLo) = (1 + Tr T) R ia (5 ab - T ab ) (R~ 1 Ry) b Jg°(\y\), (E.15) 



8tt 2 

where we have used (3.16), (A. 6), (A. 3) with W = 0, and 

tr (LtLtL*)=^He abc Jt a . 



(E.16) 



Again, (E.14) is most easily evaluated by switching to ^-integration and expressing x in terms 
of y via 

y = R- l x = U - AiR^Ry) 2 - Br 2 Tr^i?) 2 ] y - Cr 2 (R- l R) 2 y. (E.17) 

Then, the isoscalar (/-factor g I=0 of the nucleon defined by 

9i=o 



2M, 



J, 



N 



is found to be given by 

16tt eM N 



9i=o 



3V3 fn 



N 



1 + -e^p 2 (2Z + C) 
5 



(E.18) 



(E.19) 



The isoscalar mean square magnetic radius (r 2 ) MI=0 , which is defined as (E.14) with extra 
weight x 2 divided by (E.14) itself, is given by 



f °°dpp°jr(p) 


\ + le 4 n 2 N p 2 (4Z + C)] 




l + fe^p 2 (2Z + C)] 



fV (r 2 ) 

"it V /M,I=0 



Isovector magnetic moment and magnetic charge radius 
Calculation of the isovector magnetic moment, 

J d 3 x x x J 



(E.20) 



V,a=3 



>,*)» 



(E.21) 



is much more involved. Relaxing the a = 3 restriction in (E.21) to a generic a, the i-th compo- 
nent (ni=i)i is found to have three kinds of dependence on (i, a); R ia , (i?(i? _1 R) 2 ) and cuifl a . 
When we consider only the nucleon matrix elements, we have the following representation [7]: 

1 



1 1 

Ji = -jOi ® 1, 4 = 1® -T a , R ia = --CTi (X> T a . 



(E.22) 
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Then, owing to (E.13), we have 



and consequently 



-2 1 1 



{r{r- x r) 2 ). = - Wi n a - ft 2 ^ = o, 



(E.23) 



(E.24) 



both for the nucleon states. Using these facts, the isovector ^-factor g I=1 of the nucleon defined 

(E.25) 



,J 2M N 2 



is calculated to be given by 

with 
9i=i = 



9i=i 



8% M N 



9i=i, 



(E.26) 



°°dpp 2 sin 2 F<[ 1 + 4(F') 2 + 4^4^ 



+ ^ 4 ft> 2 



sin 2 F 



pZ' + 5Z - C - A{F' f (pZ' + Z + C)+ 4 — (pZ' + 3Z + 2C) 



(E.27) 



Similarly, the isovector mean square magnetic radius (r 2 ) MI=1 is given by 

, sin 2 F 



9i=i Jo [_ 



+ ^ 



sin 2 F 



pZ' + 7Z - C + 4(F') 2 (-pZ' + Z - C) + 4 — — (pZ' + 5Z + 2C) 



(E.28) 



Comparing (E.27) for gj = i with (2.18) for I and (4.2) for we find the following simple 
relationship: 



9i=i = ^ (X + e'S^j) = 



3v ^lft r 1 



(E.29) 



where the last equality is due to (E.13). This implies that gi=\ has a much simpler expression: 

2 M N 



9i=i 



ft 



(E.30) 



and that the magnetic radius (r 2 ) M J=1 is in fact equal to the electric one {r 2 ) I=l (E.ll): 

<r 2 ) M , fcl =(r 2 > fcl . (E.31) 
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The relations (E.30) and (E.31) are also valid in the rigid body approximation and should 
have a simple origin. 

Axial vector coupling 

Let us consider the axial vector current J Aa (2.8). The axial vector coupling constant 
9a = 5u(0) is obtained by identifying the nucleon matrix elements of 



o9x J A (x, t), 



with the q = p' — p — > limit of 



(N'(p')\ J\M \N(p)) = ^ 9A(q 2 ) (N'\ a, ® r a \N) 



(E.32) 



(E.33) 



valid for nonzero m n . The evaluation of (E.32) is also very complicated. Using the relation 
(E.24), we obtain 

= 53(52' + a?). (E.34) 



9a 



where g^ A and g A \ which are contributions from the part of J A linear and cubic in or 
R^, respectively, are given by 



"(3) 



5.? = 



gf = 4 / dpp sin2F<( -(F'f - F 



dp p{pF' + sin 2F + -e A h 2 N p 2 



2pF' + sin2F 4 + p 



d_ 

dp 



(2Z + C)\, (E.35) 



n/ tanF sin 2 F 



1 

5 ( 



2p 2 (F') 2 (-pZ' + 5A- 10B - AC) 



+ pF' tanF (pC + 2A + 2C + 1) + sin 2 F ^4 - (2 + p— ) (2Z + C) 
Isoscalar quadrupole moment 



(E.36) 



The isoscalar quadrupole moment operator Q\- °, 



(E.37) 



is a measure of the non-spherical deformation due to the spinning motion. Q\, vanishes 



identically in the rigid body approximation where the baryon number density is spherically 
symmetric [36]. In the present case with relativistic correction, it becomes non-trivial: 



fyQir = 



8ne' 
~15 



f°° ft 2 /1 1 \ 

J dpp" (2A - 5C) Jr(p)jt J >} ~ 3^ J2 J • (E.38) 

However, the nucleon expectation value of the operator \ { Jj, Jj} — \b~ijJ 2 and hence of (E.38) 
are equal to zero. 
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